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Abstract 

We construct the equation of a state of the classical QGP valid for all values of F = V/K, the 
ratio of the mean Coulomb to kinetic energy. By enforcing the Gibbs relations, we derive the 
pertinent pressure and entropy densities for all T. For the case of an SU(2) classical gluonic plasma 
our results compare well with lattice simulations. We show that the strongly coupled component 
of the classical QGP contributes significantly to the bulk thermodynamics across T^. 
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I. INTRODUCTION 



Classical Plasmas are statistical systems with constituents that are locally charged but 
globally neutral. An example is the electromagnetic one component plasma (OCP) also 
referred to as Jelium. A number of many body theories have been devised to analyze the 
OCP in the regime of small T = V / K , the ratio of the average Coulomb energy to kinetic 
energy l| . Most of the extensions to larger values of P are based on higher order transport 

u n 

equations [2| or classical molecular dynamics [3|]. 

The Classical Quark Gluon Plasma (cQGP) as developed by Gelman, Shuryak and Za- 
hed can be regarded as an extension of the OCP plasma to many components with non- 
Abelian color charges ^. Stability against core collapse is enforced classically through a 
phenomenological core potential. The origin of the core is quantum mechanical. Detailed 
molecular dynamics simulations of the cQGP [4] have shown a strongly coupled plasma for 
T = V/K \ oi larger. The cQGP maybe in a liquid state at moderate values of P. In a 
recent analysis [sl we have used analytical methods of classical liquids to construct the free 
energy for small P both in the dilute case and at high temperature after resummation of the 
screening effects. 

In this paper we combine the results in j4| obtained from molecular dynamics with the one- 
loop analytical results in 51] to construct the equation of state of the cQGP for all values of 
P. We will show that the strongly coupled component of the cQGP contributes significantly 
to the thermodynamics across the transition whether in the energy density, pressure or 
entropy density. In section 2, we derive the excess energy of the cQGP for small P in the 
one-loop approximation. In section 3, we use ideas from classical electromagnetic plasmas 
to interpolate between the one-loop result at low P and the molecular dynamics results at 
large P for the SU(2) cQGP. Particular attention will be given to the core parameter using 
Debye-Hiickel plus Hole (DHH) theory. In the quantum QCD plasma, P runs with T. In 
section 4, we use the interpolated excess energy density together with the Gibbs relations to 
derive the pressure and entropy of the cQGP. In section 5, we compare the results for the 
SU(2) cQGP with SU(2) lattice simulations. Our conclusions are in section 6. In Appendix 
A we show that the interaction corrections to the concentration do not affect the one-loop 
result to order pt with the bare particle concentration. In Appendix B we summarize the 
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Debye-Hiickel plus Hole theory to assess the range of the core in terms of the Debye length. 



II. EXCESS ENERGY: ONE-LOOP 



In classical plasmas the key expansion parameter at zero chemical potential isT = V/K 
the ratio of the mean Coulomb to kinetic energy. For an Abelian or QED plasma, 



while for a non- Abelian or QCD plasma [4| 



a 



r = (II.2) 

47r Taws 

with kB = I and aws the Wigner-Seitz radius satisfying N/V{4:Tcay^rg/3) = 1. C2 is the 
quadratic Casimir (C2 = q2/{N^ — 1)) and g is the strength of the coupling. In the cQGP g 
is fixed, while in the QGP g runs with temperature. The running is quantum mechanical and 
beyond the present classical analysis. In sections 4,5 it will be addressed phenomenologically. 

Since the Wigner-Seitz radius aws is tied with the density or concentration (the bare 
concentration is Cq = N/V), it is straightforward to express the free energy in terms of F. 
After resumming the screening effects and to one-loop, the free energy reads [sl 

f3 = -c ^(^c - l)l^c-^ + vr(iV^ - l)c 7 a 

-2Ti'i{Nl - lyi^ia^ + c(3^l, + 0{I3^) (II.3) 

with 7 = g'^/ATT(3C2, c the concentration and a the core radius. The concentration is deter- 
mined by the chemical equation [5] and to leading order is cq = N/V as detailed in Appendix 
A. The core a is a parameter in the cQGP much like in normal classical liquids. Its origin is 
quantum mechanical. In Appendix B, we use the classical Debye-Hiickel plus Hole (DHH) 
theory to assess the size of the core in terms of the Debye screening length. 

To see how the expansion in the concentration c in (lII.Sp converts to an expansion in F, 
we note that the Debye-Hiickel contributions (first two terms) can be rewritten as 

= -c^-ci^^(^)i(Ar,^-l)(/?C,)i 
NT N N 3 Mtt^ ^ " 
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V3 ^ o ^ 47r 0^/5 

1 , - -o . , _ 3 



-1 - ^(A^c - i)r§ (11.4) 

By defining the excess free energy F^.^ as Fe^(r) = -F(r) — -F(O) we obtain to one- loop, 



= -^{Nl - i)ri + \5{Nl - i)r^ - 3y35^(iv,^ - i)rt + o(r^) (11.5) 

with 5 = a /aws- -^(0) will be identified with the free gas or Stephan-Boltzman contribution. 
The excess energy U^x of the cQGP follows from the excess free energy F^x as 



NT Jo NT r 

For instance, the Debye-Hiickel contribution in fill. 71) yields through flII.6P the excess energy 
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-— (iv;-i)r2 (II.7) 

NT 2^ ' I J 

in agreement with the Debye-Hiickel excess energy for the cQGP initially discussed in \&\ 
using different methods. In general, the energy density splits into the free plus excess the 
contribution 

.(r).^.^ + Mn.,.,r) (II.8) 

with the free contribution eQ = ess identified with the Stefan-Boltzmann energy density e^^. 
In relative notations, 

= 1 + (11.9) 

^SB ^SB V 

Using (111.61) together with ( 111.51) . we obtain the one-loop excess energy density 



^ = -^{Nl - i)ri + \5{Nl - i)r^ - ^^V-35\Nl - i)ri + ^(r^) (ii.io) 

which is valid for small F. In Appendix A, we show that although the concentration c in F 
is not Co because of interactions, to order F^/^ we may set c = cq. 



III. EXCESS ENERGY: FULL 



In |5[ the one-loop expansion was shown to converge up to F ^ 1 for the free energy. 



see below). Larger values of F 
J]. For an SU(2) plasma (say a 



The range is even smaller for the energy with F ^ 1/2 
have been covered by molecular dynamics simulations in 
constituent gluonic plasma) the numerical results for the excess energy were found to follow 
the parametric form ^ 



mol 



NT 



1 2 2 
-4.9 - 2F + 3.2F2 + . 

F3 



(IILl) 



For Nc = 2 the one-loop result fill. 101) reads 

Uanal 



2 2 
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(III.2) 



NT 2 2 2 

To construct the full excess energy valid for all F we will proceed phenomenologically by 
seeking an interpolating formulae between (1III.1I) and (1III.2I) borrowing from ideas in classical 
plasma physics [7]. A similar approach was also advocated in Isj] using different limits. 

In the Abelian or QED plasma, the excess energy based on Debye-Hiickel theory is eval- 
uated for F < 0.1. Molecular dynamics simulations are generated for 1 < F < 180. The two 
are combined numerically through a power function in the form 7| 



'lu ,(r) 
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FIG. 1: Excess energy for QED (left) and Nc = 2 QCD (right) 
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^r<o.i(r) + /(r)Mr>i(r) 

Ufuiiin = i + f{T) ^ 

with /(r) a fitting power function of the type ar^(= 3.0 x lO^F^ ''). fllll.Sp interpolates 
smoothly between the exact analytical results at low F and the simulations at large F as 
shown in Fig. [1] (left). In the insert we show the nature of the size of the gap in the range 
0.1 < F < 1 for the Abelian plasma. 

In Fig. [T] (right) we show our = 2 results at low values of F (one-loop) and large values 
of F (simulations). The one- loop results depend on the size of the hard core a. Recall that 
the simulations in ^] are carried with a fixed higher power law repulsion to mock up the 
core. So the simulations seem to favor a small hard core in the Wigner-Seitz units. In fact, 
cr can be set by the Debye radius in the DHH theory joj] detailed in Appendix B. It changes 
with F. Specifically, 




FIG. 2: Core parameter 6 = a/aws- See text. 



s = ^ = ^^((1 + (3F)i)^ - 1) = 7:^((i + (3r)^)^ - 1) (in.4) 

aws CLwsf^D ^ ' (3F)2 ^ ' 

which is shown in Fig. [2l The core size 5 varies in the range 0.2 — 0.5 for F in the range 
0.1 — 1. We fix 5 = 0.4 in the range 0.1 — 1. With this in mind and following the Abelian 
plasma construction, we find the excess energies shown in Fig. [1] (right) to be fit by 
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UanaliV) + 5.5 X WV'-^UmoliT) .^^^ . 

^sumiT) = 1 + 5.5 ^ io^r5.4 

The power function is numerically adopted to yield a small deviation (less than 0.1%) for 
r < 0.1 and r > 1. The precise choice of the core parameter is actually not very important, 
as small changes in core size can be compensated by small changes in the power function for 
the same overall accuracy. 

IV. THERMODYNAMICS 

Knowledge of the energy density for all values of F can be used to extract all extensive 
thermodynamical quantities in the cQGP with the help of the Gibbs relations. Indeed, the 
pressure and entropy follow from the Gibbs relations 



dP 

. = ^ = . (iv.i) 

V VdT ^ ^ 

In so far the classical plasma parameter F as defined fill. 21) is a fixed parameter. However, 
in QCD it runs through a^. It is only a function of temperature. Specifically, 

r = ^ = {fco)hc2a,iT) = {0.2UiN! - l)fyC,a,iT) (IV.2) 

where 0.244(Ai"^ — 1) is the black-body concentration for adjoint gluons. The exact running 
of as{T) will be fixed below. 

Using fllV.ll) together with flIV.2p yield the pressure and the entropy density directly in 
terms of the energy density 



P 1 /"^^ f 

3_ / dT'T'^ — (T') 



PsB Jtc tsB 

s 3 e 3 1 



- — (T) + -,T^I dT'T'^ — {T') . (IV.3) 
ssB ^esB 4T^Jn esB 

Here Tc is identified with the SU(2) transition with Pc = 0. For a constituent gluonic plasma 
Tc = 215 MeV. All bulk thermodynamics is tied to the energy density by the Gibbs relations. 



V. SU(2) LATTICE COMPARISON 



To proceed further we need to know how as{T) runs with T in pure YM and QCD, to 
determine the behavior of the extensive thermodynamical quantities. The loop expansion 
allows a specific determination of the running as{T) that is unfortunately valid at high 
temperature or weak coupling. How as{T) runs at strong coupling is unknown. Here we 
suggest to extract as{T) from the lattice data by fitting our energy density ( 111.91) valid for all 
couplings to the SU(2) lattice data in [l^. In Fig. [3] we show the fit of the normalized energy 
density in the cQGP to the SU(2) lattice data. The lattice results are used with = 5 
MeV and = 215 MeV as suggested in 0. 




FIG. 3: Energy density fit to SU(2) lattice data. See text. 



The energy density fit in return yields through (]II.9|) . (jlll.SP and flIV.2|) a specific running 



of the strong coupling constant as{T) which we show in Fig. H] (left). Its corresponding 
running plasma constant F(T) is shown in Fig. H] (right). In Fig . H] we also show two running 



coupling constants extracted from lattice measurements in [HI] for comparison. Our energy 
density fit suggests lower values of the running coupling constant. 

Since we have fit the energy density to the lattice energy density, to extract as(T), it 
follows from the Gibbs relations that all the extensive thermodynamical quantities are fixed 
for the SU(2) cQGP. In Fig. [5] (left) we show the behavior of the energy density, pressure 
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T/T„ 



FIG. 4: as{T) (left) and r(r) (right) for Nc = 2. See text. 

and entropy density across Tc. In Fig. [5] (right) the trace of the energy momentum tensor is 
shown versus SU(2) lattice data. 
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FIG. 5: Bulk Thermodynamics from the cQGP versus SU(2) lattice. See text. 



Our current analysis of the bulk thermodynamics of the cQGP allows us through the 
excess energy ( IIII.5|) and the Gibbs relations (lIV.ip to assess the role of the strongly coupled 
component of the cQGP. In Fig. [6] we show the two contributions (loop and molecular) to 
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the bulk thermodynamics following from the separation 





FIG. 6: Relative contributions in the cQGP bulk thermodynamics. See text. 



/■p\ UanaliX) 
UsU{2),anal[>' ) 



1 + 5.5 X 10215-4 
5.5 X lO^T^'-^UmoiiT) 



in the energy density. The strongly coupled component of the cQGP generated by the 
molecular dynamics simulations contribute significantly across the transition temperature, 
say in the range (1 — 2.5) Tc. 
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VI. CONCLUSION 



We have constructed the energy density of the cQGP vaUd for all values of the plasma 
parameter F, that interpolates between the one- loop result at small F and the molecular 
dynamics simulations at large F. We have used it in conjunction with the Gibbs relations to 
derive the Pressure and entropy of the cQGP. 

In quantum QCD F runs through the QCD coupling constant at weak coupling. The 
running at strong coupling is unknown in general except for some recent lattice simulations 



We have suggested that a fit of our energy density to the lattice energy density [10 1 



allows an extraction o: 
lattice simulations 



Hi 



the running coupling that is smaller than the one suggested by direct 



We have used the extracted running coupling constant to predict the entropy density, 
pressure and energy-momentum trace of the cQGP. The latter compares well to direct lattice 
SU(2) simulations. We have shown that the strongly coupled component of the cQGP 
contributes significantly to the bulk thermodynamics across the transition temperature. We 
expect transport properties such as diffusion and viscosity, as well as energy loss to be also 



significantly affected in this transition region in the cQGP as we discuss next 12 |. 
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APPENDIX A: CONCENTRATION 

The bare concentration cq = N/V which is identified with the black-body radiation in 
the cQGP is in general modified to c = cq + Ac due to interactions. As a result, the plasma 
parameter F (1IV.2P is in principle different from the one used in the text. The corrected 
plasma constant is 

r. = {^cf pC^asiT) = {^{co + Ac))hc2as{T) = (yCo)^l + ^)^/3C2a.(T) 
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- (^co)^/3C2a.(T) + (^co)^^— /3C2a.(T) + 0((^)^) (A.l) 

^ 6 ^3^3 Co Co 

From [sl, the shift in the concentration reads 

c = Co + Ac = Co + cln^N^ - l){pC2)^al (T) + Oi(3^) (A.2) 
The corrected plasma constant becomes 



, 47r\ - \f?\ / 47r\ -■- s s 

= (0.244(iV2 _ ^C2a,(T) + \{Nl - 1) (o.244(Ar,2 _ \ a ^^|c,|(t) 



(A.3) 



Inserting ( ]A.3I) in the excess energy density yields 



ess 2 2 2 

- - i)(r + ^(iv.^ - i)rt)^ + \m - i)(r + ^(at,^ - i)rt)^ 

-f^m^-i)(r + ^(Ar,2-i)ri)^ 

^ -^(AT,^ - i)ri + \^{Nl - i)r^ - ^^^mK - i)rt + o^) (a.4) 

which shows that to order Pi the replacement Fc = F is allowed. 



APPENDIX B: DEBYE-HUCKEL PLUS HOLE (DHH) THEORY 

At strong coupling the Debye-Hiickel (DH) theory which is essentially a classical screening 
theory fails. Debye-Hiickel plus Hole (DHH) theory is a way to address DH shortcomings at 
strong coupling by building a hole around each charge to account for the non-penetrability 
or core in classical liquids [9] at higher density or larger F. As a result, in DHH theory of 
the cQGP a color charge density around a test charge is 
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p^{r) = { ^ ' (B.l) 



a is the size of the hole, a is a classical color index (1, .., iV^ — 1), /3 = 1/T and is the 
Debye momentum 

2 A^c-l 

The negative sign in fIB.ip reflects on the screening, with the Debye cloud left unchanged 
outside a. The hole size a is flxed by demanding that each test particle is completely screened 
through 

rfr47rrV(^) = F=Q'' (B.3) 

V47r 

This condition, flxes the size of the hole 

= -((1 + 1^)^-1) (B.4) 



In terms of 



r = ^-92_ (B.5) 



47r Taws 



the hole radius is 



a = — ((l + (3r)i)i-l) (B.6) 

after flxing the Wigner-Seitz radius a-^s through co(47ra^5./3) = 1. Again we have set c = Cq. 
From (]B.6|) it follows that the hole size is smaller the higher the density or temperature. 
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